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ABSTRACT 

The Boussinesq approximation is extended so as to explicitly account for 
the transfer of fluid energy through viscous action into thermal energy. 

Ideal and dissipative integral invariants are discussed, in addition to the 
general equations for thermal-fluid motion. 


/ 


i 



INTRODUCTION 


In this paper the equations of motion for a thermo-gravitationally 
stratified fluid are developed. These equations are pertinent to atmospheric 
and oceanic flows in which temperature variations affect fluid flow through 
boundary forces and in which fluid motion affects temperature variations 
through viscous dissipation (in addition to convection). The equations 
developed here are an extension of the Boussinesq approximation in that they 
explicitly account for the transfer of fluid energy through viscous action 
into thermal energy. These equations thus conserve physical energy, while the 
Boussinesq equations do not. 

This is not to say there are no conserved quantities in the Boussinesq 
approximation. In 2-D flows there are two: the first, sometimes called the 
energy, involves the velocity squared and the temperature (or density) 
perturbation squared. This is not strictly the physical energy, which is 
linear in the temperature perturbation. The second conserved quantity in 2-D 
flow, apparently unrecognized, is essentially the spatial correlation of the 
vorticity and temperature perturbation. 

Here the general equations for thermal-fluid motion will be discussed 
first. Then the two-dimensional version of these will be detailed. Next, the 
standard Boussinesq approximation will be discussed; in particular, the two- 
dimensional form. Finally a summary and conclusion will be given. 

GENERAL EQUATIONS 3 

In this work, the fluid under consideration will be treated as 
incompressible, i.e. , density change with respect to pressure is negligible; 
however, density will change with temperature. The fluid density can then be 


written as 



( 1 ) 


p = P 0 + p 

where p Q is a constant and p' is the density variation. Similarly, the 
temperature is given by 


T = T q + T' (2) 

where T^ is a constant and T' is the temperature variation. The relation 
between p* and T* is linear and is given by 

P' - -p 0 Pt" (3) 

where B is the thermal-expansion coefficient of the fluid (at T^). 

The pressure, in turn, is given by 

P - P Q + P^ (4) 

where p" is the pressure variation and 

P 0 = P 0 g • r (5) 

Here g is the acceleration due to gravity and r is a position vector. 

In the present case, the Navier-Stokes equation must include the effects of 
gravity : 


3v > 1 _ JL + + 

+ v. W = - “ to + V V v + g 
P v 
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( 6 ) 


(here v is the kinematic viscosity). If we now place (1), (3), (4) and (5) 

into (6) and neglect second and higher order quantities, we arrive at 

* 


3v x + * 

a* v • * ■ 


V(pVPq) + v^v - ffT'g 


This, along with the incompressibility condition 
V . v = 0 

form the velocity equations of the Boussinesq approximation^- 

Now let us consider the temperature (and by eq. 3, the density) 
evolution. (It is here that we differ from the standard Boussinesq 
approximation and use a more general form). The general equation for 
temperature evolution in an incompressible fluid is 


(7) 


(8) 
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(9) 
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where x * ~ZT~ is the thermometric conductivity (cm /sec), ic is the thermal 
P 

conductivity, and is the specific heat (at constant pressure) of the 
fluid. The last term in (9) is 


i jv 2 Jii iii , \ 
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where the summation convention is in effect (e.q., + v 2 w 2 + 

v 3 w 3 ^ * 


At this point, let us write the continuity equation: 

-|£ + 7 . pv = 0 

Expanding the divergence term yields 
+ v . 7p - -p7 . v 

Using (1) and keeping only the lowest order terms (remember 

that 7 . v *0) gives 
3p" -► _ ♦ 

ir + v • 7p - -p 0 7 • v 

Now, because of (2), equation (9) may be written 


(ID 


(12) 


(13) 
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(14) 


Using (3), a comparison of (13) and (14) yields 


-> 3v. 3v. 2 

v . 8 [xvV . ~ i <15) 

The term within brackets in (15) and on the right hand side of (14) is small 

compared to the convective term (v . VT ') ; also, the thermal expansion 
-2 


coefficient (K10 


The product of $ and the bracketed term in (15) is 
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therefore of a higher order than any terms being retained in this theoretical 


development. It is for this reason that we accept the incompressibility 
condition (V • v = 0), i.e., any possible contribution of V % v to equation 
(7) is negligible. These assertions will be made quantitative in the next 


section. 

Let us now examine the energy conservation implicit in (7), (8) and 
(14). To begin, we use the vector identity 

vx ( Vxv) = — W - v . W (16) 

to change (7) into 

— = v x a) - 7(p7p 0 + V /2) + vv v - ?T^g (W) 

where the vorticity is to = V x v . Now, (17) is multiplied by v • , (14) by 
and the resulting equations are added, resulting in 


-!- (- v 2 + C T ') = -V . [v(^ + + C T') - C 

St v 2 p P Q 2 p p* 

+ v[v . + ( ^ v i )(3 i v k )] 

( 1 € 

- [ ffT v . g] 

(Here S/Sx^ = 9^. for convenience). At this point we must treat the last two 
bracketed terms on the right hand side of (18) in more detail. 


First, consider the term whose coefficient is v : 


V . fv + 


1 W W + (3 kV (3 iV 
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Here V . v = 3^ = 0 is the incompressibility condition, 

the last term on the right hand side of (18) is 
- 0T‘ 


Using g 


= - V* 


-► 

V 


+ p' -► + 

g = — v . g 


o + + + 

o v • v < Pn9 • r) 


T v • % 


pp 0 - 

v . ( — - v) 

p o 


— v . ( p"v) 


( 20 ) 


Using the continuity condition, (20) becomes 


- 0T" v 


g = V 


. P P 0 +. P 0 3p" 

. ( — r— v) + — — — K — 

p 0 p 0 * 
pp, 


It (-2j PT * V . <-y 


v) 


- - -T 2 [4r (P n T") + v . (p T' v)] 
Pq 3t 0 0 

Now, using (19) and (21), equation (18) becomes 


( 21 ) 


-v . [v (^ + T- + {C + p } T "*) - — VT' - v(-^Vv 2 +v . Vv)] 


p o‘° 


( 22 ) 


This energy conservation equation is accurate up to first order in terms such 

as p" and T Multiplying both sides of (22) by p Q and integrating this 

expression over the volume of interest, V, we arrive at the following 

— E = - ff F . dS (23) 

dt JJ S 


where the total energy in volume V is 

E - I, [ T “o + ">o c p * T ' ] *' 


(24) 
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and the energy flux is 
F = v [p' + p„ V 2 /2 + (C + 5 pj T '] 

0 p 0 

- ic vr' - n (“ Vv 2 + v . V v) (25) 

where the viscosity is n = p Q v. Note that, in (24), 

1 2 

— p^ v : Kinetic energy density 

p C T ■*: Thermal energy density 

0 p 

8 PgT * = -p^gz: Gravitational potential energy density (26) 

Several points should be noted concerning the energy conservation 
formulas (23), (24) and (25). First, the reference point for the total energy 
is arbitrary; we can add a constant Eg to E in (24) without changing (23). 
Second, if the volume V is large enough so that there is no energy flow 
through the bounding surface (or if the flow is periodic or balanced across 
the boundaries) then the surface integral in (23) is zero and dE/dt = 0. 

In the next section, the equations will be written in non-dimensional 
form. This will allow for the relative size of various terms to be determined 
and establish criteria for deciding on the similarity of flow in different 
physical situations. Application of these results to flows within the 
NASA/LaRC Vortex Research Facility will provide a quantitative example. 

MONDIMENS IONAL EQUATIONS 

At this point we will introduce the presence of a temperature gradient. 
The fluid, at rest, will have a uniform vertical temperature gradient 
Y (°C/cm); the temperature variation will thus be composed of two parts: 

T "* = YZ + T (27) 

Here t is the temperature anomaly field. Placing (27) into (14) gives 
4?+ v . Vt=-v.(yz) + X^ 2 t + -r~ ( 3. v + 3 v ) 2 (28) 

at 2C K x 1 K 

P 
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This, along with the equation for fluid motion (7) and incompressibility 
condition (8) are the equations of motion for our thermal-fluid system; these 
are the equations we wish to put in nondimens ional form. 

In order to produce our nondimensional equations, we will use the 
following characteristic values: 

Table I Characteristic Values 


L: length 

V: velocity 

ys temperature gradient 

B: coefficient of thermal compressibility 

g: gravitational acceleration 

v: kinematic viscosity 

X: thermometric conductivty 

C p: specific heat at constant pressure 

•« *• •• 

N = /Bgy: Brunt-Vaisala frequency 

4 

The quantities in Table I can be combined into nondimensional groups ? 
these groups have values for specific flows which allow them to be compared to 
other flows. The pertinent groups are given in Table II. 
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Table II Dimensionless Groups 


Reynolds number: 
Prandtl number: 
Strouhal number: 
Froude number: 
Grashof number: 
DuLong number: 


P = — 


VL 

v 

v 

x 

V 

NL 

V£ 

Lg 


1 

N* 


Du = 


( RN* ) 

V 2 
VL C 


We may now take equations (7) and (14) and transform them into 


dimensionless equations. In these equations, v will be measured in terms of 
V, T^ in terms of YL, (pVp^) in terms of gL, and time t in terms of N \ 


Remembering this, the equations are 


1 ■*• + 1 

9t + N* V W “ FN* 


V (*-) + 


1 

RN* 


V 2 ^ 


N* T " z 


(29) 


9T 1 + 1 

3t + N* v * N*RP 




(30) 


In order to estimate the values of the coefficients in (29) and (30) for 
a practical situation, we have drawn the characteristic values from a typical 
experiment at the Vortex Research Facility at NASA/LaRC^* These are presented 
in Table III, along with the associated dimensionless numbers. 

At this point we can estimate the various sizes of the terms in (29) and 
(30). The right hand sides of (29) and (30) are important in relation to the 
size of the v • V terms on the left sides, which have coefficients of 1/N*. 
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Table III Typical Values 


V * 50 cm/sec 

L *» 10 cm 

v~ .145 cm 2 /sec 

X * .202 on /sec 

y = .01°C/cm 

g * 980 cm/sec 2 

$ = 3.48 x 10- 3 /°C 

C p - 1.012 x 10 7 cm 2 /(sec 2 ”°C) 

N = 0.185 sec -1 
R = 3450 
P = 0.72 
N* = 0.037 
Du =* 0.0025 
F * 0.255 
G = 16300 

Thus, the relative sizes of the various terms to these "convective” or 
inertial terms is given by their coefficients with the 1/N* factor divided 
out. These relative coefficients are given in Table IV. 

Table IV Relative Coefficients 

v equation 

i * 3- 92 

F 

— * 2.9 xio" 4 

R 

(N*) 2 = 1.4 x 10 


T ^equation 

— = 4.0 x 10" 4 
RP 


7.2 x 10" 7 
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At this point, let us estimate the validity of assuming 


incompressibility. If we nondimensionalize (15) we have 


7 • * ‘ ,ar ' + 2 ( 8 k V i + S i V k> 2 

P 


( 1 .4 x 10“ ? ) V^T " + (2. 


(2.5 x 10 8 ) ( 3 k v ± + v^) 2 


Thus, in a dimensionless representation where all variables and their 
spatial derivatives are nominally of unit magnitude, the error in assuming 
incompressibility is about one part in ten million. 


TWO DIMENSION AL FLOW 

It is often useful (and computationally expedient) to approximate certain 
flows as being two dimensional. In particular, the evolution of vortex wakes 
trailing aircraft and submarines can be studied by examining the dynamics of 
the pertinent fluid in a cross-sectional plane transverse to the direction of 
motion. In this section we will reduce the general equations of motion (8), 
(17) and (28) to those appropriate to two-dimensional thermally-stratified 
flow. 


As is well known, the incompressibility condition (8) can be 

automatically satisfied in 2-D flow by expressing the velocity in terms of a 

stream function ip (y,z): 

* 

v Vx(xip) 


= - x x V ip 

_ ; m . : it 

' y 3z Z 


(31) 


The vorticity, in turn, can also be expressed in terms of the stream function: 



( 32 ) 


u =Vxv=Vx(VxxiJ>) 

A A a 

= V(x . Vi|>) - xv 

- - X T 2 

The scalar vortlcity is thus w = -V 2 where the operator 
2 2 2 2 2 

7^ = 3 / 3y + 3 / 3z because there is no x dependence in the 2-D 


approximation . 

At this point let us take the curl of eq (17): 

= Vx(v x u) + vV 2 to - B7T' x g 

Similarly, the divergence of (17) gives 

^ (^— + v 2 /2) = 7. (vxu)-BVr'.g 
P 0 


(33) 


(34) 


These equations, (33) and (34), are still applicable to 3-D 
incompressible flow; in particular, they show that the pressure p' is given, 
at each instant of time, by the instantaneous values of v and T ' , both of 
which obey non-linear coupled evolution equations. If we now put (27), (31) 
and (32) into (28), (33) and (34), we get a set of equations for studying the 
2-D flow of a thermally-stratified fluid: 


-If- 


ll = ( J4 Jl _ Ji Jl) + VV 2 T + Y -P 

9t ' 3y 3z 3z 3y* X T 3y 


a , 2 a , 2 


^(p'/p 0 ) “ hj 2 - •“ ^,v 2 , + Bg( Y + 


( 35 ) 


(36) 


(37) 


We also note that the quantity conserved by these equations, for a closed 
system, comes from (24); this is the energy of the system (leaving off time 
independent terms ) : 
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E 


(38) 


* / v2 + (p o c P + 8 p o )T]dydz 

V 

We will call this assembly of equations, (35), (36) and (37), the thermal- 
vortex equations* 

Let us note here that the energy given by (38) is conserved in the 
presence of viscosity and thermal conductivity. As has already been pointed 
out, the energy as given by (38) is explicitly composed of a kinetic part, a 
thermal part, and a gravitational part. This conserved quantity (integral 
invariant) E can thus be considered to be the actual physical energy of the 
fluid system we are modelling. 

In comparison, let us look at the system of equations which are obtained 
when we set v * x 38 0s the inviscid equations. It is well known that in this 
case, (35) and (36) conserve a certain quantity which has been denoted the 
•energy* 2 ; this energy E^ is 

E' = J v P Q (^v 2 + T 2 )dydz (39) 

This, however, is not the physical energy (although, in some computer 
models, it will be numerically close). In addition, with v * X s * 0 / it is 
easy to show that, for a closed system, there is also another conserved 
quantity: 

C = / v tax dydz (40) 

The presence of two conserved quantities fundamentally changes the 
evolutionary behavior of a dynamic system 6 . If the physical system which is 
being modelled does not actually have two (or more) integral invariants, 
although our model system does, then the predictions we make with our model 
system must surely be suspect. 
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These model systems arise in a number of applications, notably inviscid 
hydrodynamics and magnetohydrodynamics. These are interesting models; their 
numerical realizations indeed give rise to very interesting statistical 
systems whose fundamental behavior depends on the number of integral 
invariants 7 . However, viscosity (and thermal conductivity) are always present 
in real systems so that the accurate modelling of any real system requires the 
incorporation of the proper dissipative mechanisms (viscosity and 
conductivity) . 

The presence of integral invariants allows us a numerical check of our 
simulations: we set v ■ x = 0 and see that the integral invariants are 

conserved: The real system has v * 0, x * 0; in this case there is only one 

"integral invariant": the total energy. 

COHCLOSIOH 

In this paper, the equations of motion suitable for studying the 
incompressible motion of a thermally stratified fluid were developed. It was 
shown that the two quantities which satisfy time evolution equations were the 
temperature anomaly and the vorticity. The equations were therefore termed 
the thermal-vortex equations . 

These expressions (for a closed system) contain one integral invariant, 
which is clearly the physical energy. In contrast, the quantity commonly 
called the energy of the non-dissipative equations (v = x = 0) is not, in 
fact, the physical energy. In addition, the non-dissipative equations contain 
a second integral invariant. The presence of two invariants, plus the fact 
that neither is the actual energy, divorces the non-dissipative model from 
physical reality. 
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